Coulomb gap and variable range hopping in a pinned Wigner crystal 
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It is shown that pinning of the electron Wigner crystal by a small concentration of charged 
impurities creates the finite density of charged localized states near the Fermi level. In the case of 
residual impurities in the spacer this density of states is related to nonlinear screening of a close 
acceptor by a Wigner crystal vacancy. On the other hand, intentional doping by a remote layer 
of donors is a source of a long range potential, which generates dislocations in Wigner crystal. 
Dislocations in turn create charged localized states near the Fermi level. In both cases Coulomb 
interaction of localized states leads to the soft Coulomb gap and ES variable range hopping at low 
enough temperatures. 



The growth of GaAs-heterostructures mobilities has 
made possible to study collective properties of the hole 
gas with a very small concentration. The average dis- 
tance between holes r s , calculated in units of the hole 
Bohr radius, has reached 80 and exceeded the theoreti- 
cally predicted freezing point of the hole liquid, r s = 38, 
where it should become the Wigner crystal. The trans- 
port of such low density systems with and without mag- 
netic field was studied recently. If a strong magnetic field 
applied along the hole plane the low temperature conduc- 
tivity of the Wigner-crystal-like phase seems to be related 
to the variable range hopping (VRH) 1 . Thus, it seems 
timely to discuss ohmic transport of the Wigner crystal. 
Of course, any Wigner crystal is always pinned by im- 
purities and can not slide as a whole in a small electric 
field. Interstitials and vacancies of the Wigner crystal 
provide conductivity but they have large activation en- 
ergies and freeze out at low temperatures. On the other 
hand, the VRH conductivity requires nonzero density of 
states of charge excitations near the Fermi level, both in 
the case of Mott law or Efros- Shklovskii (ES) law. If pin- 
ning centers produce such "bare" density of states then 
the long range Coulomb interaction creates the Coulomb 
gap around the Fermi level and leads to the ES variable 
range hopping. In this paper, I study the origin of the 
bare density of states of the pinned Wigner crystal, the 
Coulomb gap on its background and the resulting VRH 
conductivity at low temperatures. 

Let me start from a clean electron Wigner crystal on 
uniform positive background. Consider an additional 
electron added to the perfect Wigner crystal and let the 
lattice relax locally. This addition costs the interstitial 
energy 0.14e 2 /«;a, where a is the lattice constant or the 
optimal triangular lattice, k is the dielectric constant 
of GaAs^. On the other hand, the energy cost of the 
extraction of an electron and the corresponding relax- 
ation around the vacancy is equal to the vacancy energy 
0.26e 2 / na. Thus, the density of states of relaxed exci- 
tations (electronic polarons^), consists of two delta-like 
peaks at energies — 0.26e 2 /«a and 0.14e 2 / na. The lower 
peak corresponds to occupied states, the upper one con- 
tains empty states. Between them there is the hard (com- 
pletely empty) gap of the width 0.40e 2 //ca. The Fermi 
level of electronic polarons (the zero of energy) is situated 



in this gap, therefore the density of states at the Fermi 
level is zero. 

In typical samples the Wigner crystal is pinned by 
charged impurities: donor and acceptors. Even if there is 
no intentional doping and the two-dimensional electrons 
gas (2DEG) is created in a heterostructure by a distant 
metallic gate, there is always a small concentration of 
residual donors and acceptors. Many high mobility het- 
erostructures are intentionally doped by a donor layer, 
which is separated from 2DEG by the spacer with a large 
width s 3> a. This is an additional source of disorder. 
In the beginning of this paper we neglect this source and 
deal with residual impurities only. We return to the role 
of doped layer in the second part of this paper. 

Let me argue that close residual donors and acceptors 
create localized states at the Fermi level. I use theory of 
pinning of the Wigner crystal by charge impurities de- 
veloped in Ref. 0. As in that paper I talk about Wigner 
crystal of electrons and concentrate on close residual ac- 
ceptors which according to Ref. produce stronger pin- 
ning than donors. (For the hole crystal the role of accep- 
tors is of course played by donors.) Because the spacer 
width, s, and the lattice constant of the Wigner crystal, 
a, are much larger than the lattice constant of AlGaAs 
the authors of Ref. |4j assumed that acceptors are ran- 
domly distributed in space with three-dimensional con- 
centration N. Each acceptor is negative because it has 
captured one electron from 2DEG. 

It was found 4 that effect of an acceptor dramatically 
depends on its distance, d, from the plane of the Wigner 
crystal. If d is larger than 0.68a the electron crystal ad- 
justs to acceptor in such a way that in the ground state 
an interstitial position of the crystal is exactly above the 
acceptor. In this case, interaction of the acceptor with 
the Wigner crystal can be calculated in the elastic ap- 
proximation. On the other hand, when d is smaller than 
0.68a acceptor creates a vacancy in the crystal, which po- 
sitions itself right above the acceptor. This means that 
the close acceptor effectively builds in the crystal cite 
pushing away to infinity its electron. We can consider 
the latter acceptor as an empty localized state for elec- 
tron, while the former acceptor can be considered as a lo- 
calized state occupied by an electron. Their energies are 
equal at d = 0.68a, so that all acceptors with d > 0.68a 
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are occupied and all acceptor with d < 0.68a are empty. 
Thus, acceptors with d = 0.68a are at the Fermi level 
and the bare density of states near the Fermi level 
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Na 
e 2 1 ' na 



(1) 



where Na is an estimate of the two-dimensional con- 
centration of acceptors located within distance a (only 
they provide substantial pinning) and e 2 / na is the char- 
acteristic energy of their pinning. The nonzero density of 
states gB makes possible VRH conductivity of the pinned 
Wigner crystal. For example, an electron from an accep- 
tor with d = 0.68a + can hop to a distant one with 
d= 0.68a -0. 

The long range Coulomb interaction between localized 
states leads to the soft Coulomb gap at the Fermi level^. 
Shape of the Coulomb gap depends on the interaction, 
U(r), of the hopping electron in the final state of the 
hop and the empty place it has left behind (excitonic 
term)2i£. If the interaction can be described by the stan- 
dard Coulomb law Uc(r) = —e /nr, where r is the two- 
dimensional distance between initial and final state, then 
the Coulomb gap has the standard for two-dimensional 
case form 
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The attraction potential U(r) generally speaking is dif- 
ferent from Uc{f) due to screening of the Coulomb in- 
teraction by elastic deformations of the Wigner crystal 
or in other words by its polarization. In order to cal- 
culate Fourier image U(q) of U(r) we have to introduce 
the Larkin length L at which acceptors destroy the long 
range order of the crystal. In the two-dimensional Wigner 
crystal pinning by close acceptors is so strong that this 
length is close to the average distance between them, 
L ~ (Na)^ 1 / 2 . Using L, we can write 



U{q) =2w/q K e(q), 



where 



qr D g z 



(3) 



(4) 



is given by interpolation between large and small q cases. 
Here rn ~ —0.32a is the linear Debye screening radius of 
the Wigner crystal. (The asymptotics of Eq. 10} at large 
and small q are similar to those of the dielectric constant 
of a two-dimensional electron gas in a strong magnetic 
fieldSi^, where the Larmour radius plays the role of L.) 

In the real space the Fourier image of Eq. Q results 
in the Coulomb potential U(r) = Uc{t) only at r > 
L 2 /a — 1/Na 2 . One can say that by these distances all 
electric lines of a probe charge located in the reference 
point in the plane of Wigner crystal leave the plane. This 
leads to the standard Coulomb gap Eq. @, but only 
at energies \E\ < e 2 Na 2 /n. At smaller distances r -C 



1/Na 2 potential U(r) is weaker than Uc{r) and grows 
only logarithmically with decreasing r. But this does 
not lead to faster (exponential) growth of the density of 
states at \E\ ~ e 2 Na 2 /n as it would if gs were very large. 
The reason is that at energy \E\ ~ e 2 Na 2 / k the Coulomb 
gap density of states, Eq. @, already reaches the bare 
density of states, gs, given by Eq. Q. The Coulomb 
gap is just a depletion of the density of states on the 
background of the bare density. Therefore, the density 
of states in the Coulomb gap can not be larger than gg£- 
This means that at \E\ ^> e 2 Na 2 / 'k the density of states 
saturates at the level of Eq. (JJ. 

The role of distance L 2 /a — 1/Na 2 can be interpreted 
in another way. When we add an electron to the clean 
Wigner crystal it makes an interstitial but its charge 
spreads to infinity. In the a pinned Wigner crystal an 
added charge is smeared in the disc of the finite radius 
R. We can find R minimizing the sum of the Coulomb 
energy e 2 / kt of the disc of radius r and the shear energy 
/j,(u/L) 2 r 2 necessary for the disc dilatation by the area 
a 2 . Here /z ~ e 2 /a 3 is the shear modulus and u ~ a 2 /r 
is the necessary displacement. It is important that the 
characteristic distance of the variation of the shear dis- 
placement (which percolates between strongly pinning 
acceptors) is the average distance between close accep- 
tors L ~ (iVa) -1 / 2 . Minimizing this sum we find that 
R = L 2 /a = 1/Na 2 . This again tells us that at distances 
larger than L 2 /a = 1/Na 2 we deal with the unscreened 
Coulomb potential. 

At low enough temperatures the Coulomb gap leads to 
Efros-Shklovskii(ES) law of the temperature dependence 
of the variable range hopping^ 



a = ( 7oCxp[-(T BS /T) 1/2 ] 



(5) 



Here <7o is a prefactor, which has only an algebraic In- 
dependence and 



T ES = Ce 2 M, 



(6) 



where C is a numerical coefficient close to 6 and £ is 
the localization length for electron tunnelling with energy 
close to the Fermi level in the Wigner crystal. We will 
discuss the value of £ in the end of the paper. 

When with increasing temperature the width of the 
band of energies used for ES hopping (TTes) 1 ^ 2 reaches 
e 2 Na 2 /k ES law is replaced by the Mott conductivity 



a = <jq exp[ 



-{T M /Tf/% 



(7) 



where Tm = Cm I '{9b£, 2 )- Transition from the ES law to 
the Mott law happens at T - T| s /T^ - T ES (Na 2 £) 2 . 
Thus, while ES law does not depend on the acceptor 
concentration N the range of ES law shrinks when con- 
centration of acceptors decreases. 

This example emphasizes universality of the Coulomb 
gap and ES law. The Coulomb gap was derived for lightly 
doped semiconductors, where disorder is as strong as in- 
teractions^. Later EfrosS suggested a model, where elec- 
trons are located on sites of a square lattice with +1/2 
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charges, the number of sites being twice larger than num- 
ber of electrons. Random energies of sites are uniformly 
distributed in the band A (they are measured in units of 
Coulomb interaction of two electrons on nearest sites). 
In this model the Coulomb gap does not survive in the 
small disorder case A << 1 because in this case in spite 
of small disorder positive charges of empty sites and 
negative charges of occupied sites alternate in the per- 
fect NaCl-like order. One could, therefore, say that the 
Coulomb gap is the property of strong disorder only. The 
above example of the Wigner crystal on the uniform back- 
ground pinned by rare strong impurities however shows 
that such impression may be misleading. Even in the case 
when concentration of acceptors N is small (one can call 
it the weak disorder case) the Coulomb gap and ES law 
survive. 

One can interpret what happens with the Wigner crys- 
tal in the terms of the Efros model. At strong disorder 
when Na 3 ~ 1 our model is close to the Efros model with 
A ~ 1. When parameter Na 3 becomes much smaller 
than unity and the disorder becomes weak, the energy 
scatter of states created by impurities stays at the level of 
e 2 /«;a, while the long Coulomb interaction between these 
states becomes much smaller, of the order of e 2 {Na) 1 / 2 
kappa. Introducing a two-dimensional random lattice 
with cites occupied by acceptors we can come to a renor- 
malized Efros model with A ~ l/(Na 3 ) 1/2 > 1. This 
brings us again to the Coulomb gap at small energies. 

Until now we talked about the role of residual accep- 
tors. Let us qualitatively discuss intentionally doped 
heterojunctions where donors are situated in the narrow 
layer parallel to the plane of the 2DEG (delta-doping) at 
a large distance s 3> a from 2DEG. Random distribu- 
tion of donors in this layer creates fluctuations of their 
potential of all wavelengths, but only harmonics of the 
random potential with wavelengths larger than s reach 
2DEG. Let us assume that the two-dimensional concen- 
tration of charged donors in the doped layer is N2- In 
principle N2 can be different from the two-dimensional 
concentration of electrons in the Wigner crystal n (2DEG 
can be compensated by acceptors or created by illumina- 
tion). If N2 <C n the Wigner crystal screens external 
potential by small, purely elastic deformations. Due to 
these deformations energies of an interstitial and a va- 
cancy depend on a coordinates. Therefore, both delta- 
like peaks of the density of states are somewhat smeared. 
Their width, however, is much smaller than the hard gap 
between them. Thus, the hard gap is preserved and no 
states at the Fermi level appear. When the concentration 
of donors, N2, grows, the amplitude u of displacements 
of electrons in the Wigner crystal with wavelength of the 
order s reaches the lattice constant of the crystal a. Such 
large deformations resolve themselves by creation of dis- 
location a 10 ! 11 , because the energy price of the dislocation 
core becomes smaller than an elastic energy, which is 
eliminated by the dislocation. This happens when the 
donor concentration, N2, is of the order of electron con- 
centration, n. Indeed, if we cover the layer of donors by 



squares with the side s the typical fluctuation of number 
of charges in a square is equal (-/V 2 s 2 ) 1//2 . Potential of the 
charge of this fluctuation reaches 2DEG practically with- 
out compensation by oppositely charged fluctuations in 
neighboring squares. As a result at n — N2 2DEG has to 
provide (Ns 2 ) 1 ^ 2 — s/a new electrons to screen random 
potential in a square. This may be done by an addi- 
tional raw of electrons in the square or, in other words, 
by two dislocations. Thus approximately one dislocation 
appears in a square with a side s. This picture is actually 
a simple visualization of the Larkin domain. Below we 
concentrate on the case N ~ n. In this case, the length 
s plays the role of the Larkin length L. 

An isolated dislocation brings the electronic polaron 
state right to the Fermi level. Indeed, if we add an elec- 
tron from the Fermi level to the end of an additional raw 
terminated by a dislocation and let the rest of electrons 
relax, the dislocation just moves along this additional 
raw by one lattice constant a and the energy remains 
unchanged. Similarly, if we extract an electron from the 
end of the raw to the Fermi level, the dislocation moves 
in the opposite direction and the energy does not change. 
This means that if we neglect the interaction of disloca- 
tion with long range ptential of donors and interaction 
between dislocations they create a delta peak of the den- 
sity of states right at the Fermi level. This peak is nor- 
malized on the concentration of dislocations. The Fermi 
level is pinned in the middle of this peak. 

In the long range fluctuating potential of donors the 
peak of the density of dislocation states is smeared, be- 
cause dislocations strongly interact with the gradient of 
potential. The peak is smeared also due to the interaction 
between the dislocations. When an additional electron is 
absorbed by a dislocation and, as a result, the disloca- 
tion moves by one lattice site, the logarithmic interaction 
with other dislocation changes. 

To understand the role of Coulomb interaction in the 
density of states of the pinned Wigner crystal one has to 
concentrate on the fate of the charge of an added electron. 
In the clean Wigner crystal, if we have a single disloca- 
tion and move it by one lattice constant adding a new 
electron, the charge of this electron spreads to infinity, 
leaving the dislocation neutral. 

If we have a gas of dislocations in the positions fixed by 
disorder and their interaction, a single charge can only 
spread to the finite distance, R. It can be estimated in 
the way we did for residual randomly distributed accep- 
tors. If we assume that dislocation are fixed in space 
by fluctuating donor potential the extra electron charge 
spreads optimizing the sum of its Coulomb energy and 
the energy of the shear deformation. Optimization leads 
to R <~ s. Thus, the Coulomb potential of charges in- 
side the pinned Wigner crystal is valid at distances in 
the plane r 3> s. 

This system is clearly similar to the Coulomb glass 
and therefore has the Coulomb gap. Indeed, as we men- 
tioned above the derivation of the Coulomb gap is based 
upon the observation that when a localized electron is 
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transferred to another localized state one should take 
into account its 1/r Coulomb attraction with the hole 
it has left (the excitonic effect). We claimed above that 
that the Coulomb interaction is valid if r 3> s. There- 
fore, the linear in energy Coulomb gap appears at the 
Fermi level. The width of this gap is e 2 / kR — e 2 /ks. At 
s 3> a the Coulomb gap occupies only small fraction of 
the energy range between interstitial and vacancy peaks. 
Away from the Coulomb gap the density of states is al- 
most constant. Thus, in the case of intentional (5-doping 
by remote donors we again arrive to the Coulomb gap of 
density of states and correspondingly to ES law at low 
temperatures. 

Let us discuss the value of localization length £ in 
Eq. ©. In classical Wigner crystal £ ~ ajr\ . This 
is a quite small value which leads to too large To and 
very large resistances in the range of ES law. However, 
close to the melting point of the Wigner crystal £ can be 



much larger, making an observation of ES law in pinned 
Wigner crystal more realistic. Magnetoresistance, ob- 
served in Ref. G] may be related to ES variable range 
hopping conductivity. 

In conclusion, I emphasize again that I am dealing with 
the Wigner crystal which in the absence of impurities 
slides on the positive background. Impurities pin this 
Wigner crystal and lead to ES VRH conductivity due to 
hops between its pieces. This situation is similar to what 
happens in a system of many quantum dots situated in 
random electrostatic potential of stray charges or in a 
system of many electron puddles with random positive 
charge of their background. Similar physics was recently 
theoretically studied in quasi-one-dimensional systems^. 
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